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14. Resumo/Notas 





Experimentos de previsao de safras em grande escala qua utili 
zam imagens de satelite, tuis c<.mo o projeto LACIE (feito nos E.U.A.), indiia 
ram a importancia de procedimentos de classificacao adequados ; tais p: 'Oce_ 
dimentos ubjetivam a medida da proporqdo de culturas agricolas em segmentos 
de imagens de satelite, com a consequents agregaqdo dos re suit ados em urn sis_ 
tema integi'ado de previsao de safras. Com este objetivo, vem sendo desenvolvi_ 
do no INPE urn pvocedimento inievalivo que envolve analista e computador, com 
quatro fases distintas: a) rotulaqdo e trdnamento ; b) classificacao; c) esti 
mativa de proporqao ; d) avaliaqao e aprovaqdo. Este procedimentc tan ainda 
as seguintes caracteristicas : 1) uso de varias passagens de satelites sobre a 
mesma area (classificacao multitempoi'al) ; 2) necessidade de urn mtnimo de in 
formaqao de campo; 3) capacidade de vevificaqdo entre classificacao automati 
ca (feita por computador) e rotulaqao do analista. Neste trabalho apresenta- 
~se uma descricdo do fluxo ih processamento, ben corns discutem-ce os prir.cipais 
algoritmos envolvidos. Dar-se-a enfase, ainda, aos resultados obtidos e as 
perspectives futuras dessa tecnologia. 
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1. INTRODUCAO 


0 Institute) de Pesquisas Espaciais (INPE) vem desenvo.lven 
do um si sterna de estimativa de safras agricolas para as culturas de ca 
na-de-acucar, soja e trigo, a nTvel nacional, baseado, principalmente, 
em informacoes colhidas por satelites. As atividades do sistema de pre 
visao de safras sao as seguintes: 

- estimativa de area plantada; 

- estimativa de prooutividade; 

- integracao da area e produtividade para a obtencao da producao 
das culturas. 

A metodologia proposta para a estimativa de area de uma 
dada cultura, para fins de previsao de safras, i a de amostragem estra 
tificada, que envolve os seguintes passos: 

1) A regiao em estudo e dividida em zonas homogineas, distribuTdas^ 
de acordo com caracterTsticas de clima, solos e tecnicas agrT 
colas comuns. 

2) Cada zona homogenea e particionada em estratos. 

3) Para cada estrato sao escolhidas, de maneira alestoria, unidades 
amostrais denominadas segmentos. 

4) A proporcao da cultura em cada segmentr 5 determinada. 

5) As proporcoes estimadas para os segmentos sao agregadasee esU 
mada a area da cultura no estrato. 

6) As areas dos estratos obtidns no passo 5 sao agregadas e assim 
e obtida a area da cultura para a zona; a uniao das areas esU 
madas para as zonas fornece a area total da cultura para a re 
giao. 
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0 presente trabalho apresenta um procedimento para deter^ 
minagao da proporgao de cultura nos segmentos (unidades amostrais). 0 
procedimento devera envoi ver a participagao de analistas (agronomos, fo 
to-interpretes) e computadores, e funcionara de modo interativo.com qiM 
tro fases distintas: 

- rotulagao e treinamento; 

- classificagao; 

- estimativa de proporgao; 

- avaliagao e aprovagao. 

Este procedimento tern as seguintes caracterTsticas: 

(a) Com o intuito de aumentar a discriminabilidace da cultura, o 
procedimento usa varias passagens do satelite. Para isto, seg 
mentos obtidos em varias datas distintas serao digitalmente so 
brepostos ("registrados"). 

0 

(b) 0 uso de informagao de campo ("verdade terrestre") devera ser 
minimizado, em fungao do alto custo geralmente associado a e£ 
te tipo de dado. 

(c) 0 procedimento dispoe ainda de capacidade de verificagao entre 
a classificagao automatica (feita por computador) earotulagao 
do analista. 

Experimentos anteriores, como o LACIE (Experimento para 
Previsao de Safras de Grandes Areas), realizado nos E.U.A., evidenciaram 
a importancia do uso de metodos adequados para estimativa de proporgao 
das culturas nas unidades amostrais. 0 procedimento desenvolvido e basea 
do na experiencia do LACIE (Heydorn et alii 1978), como tambem em traba 
lhos anteriores realizados no INPE, nas areas de Reconhecimento e Classic 
ficagao de Imagens. 
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2. DESCRICfiO DO PROCEDIMENTO 

0 Procedimento de determinagao da proporgao de cultura nos 
segmentos consiste em quatro fases distintas: rotulagao - treinamento, 
classifigagao, estimagao da proporgao de eultura no segmento e avalia 
gao-aprovagao, com o fluxo de processamento ilustrado na Figura 1. 


ROTULAQAO/ 

TREINAMEN- 

TO 


CLASSIFICA 


ESTIMATIVE 

/ DA 

AREA 



AVALIAQAO 


AGREGAQAO 


Fig. 1 - Fluxo de processamento. 


2.1 - ROTULACAO - TREINAMENTO 


Nesta fase, o analista manualmente rotula pon+os de ' ima 
gem, com dois objetivos: (1) fornecer observagoes necessarias para estji 
mar parametros para a classif icagao; e (2) dar subsTdios para testar a 
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qualidade das estimativas. Para a rotulacao, sao escolhidos aleatoria 
mente dois subconjuntos d a pontos em urn segmento (denominados tipo 1 e 
tipo 2); os pontos do tiDO 1 sao utillzados para "treinar" o classifica 
dor e os do tipo 2, para a estimativa de erro e correcao de tendencia 
da classificacao. 

0 analista podera atribuir apenas dois rotulos aos pon 
tos: cultura "x" (por exemplo, soja) e "outros" (exemplos, pastagem ); 
estas tambem serao as classes resultantes da classificacao. Alguns re 
finamentos sao introduzidos nesse procedimonto proposto pelo INPE, em 
relacao a "Procedure 1" do LACIE (Heydorn et alii, 1978). rtais especj^ 
ficamente, as imagens a serem rotuladas poderao, a criterio do analis^ 
ta, nao correspoder as cenas originalmente geradas pelo satelite 
LANDSAT, mas sim a novos gerados a partir das imagens originais por 
trcnsformacoes e filtragens; duas possibil idades se apresentam: a obten 
cao de novos canais representatives de "brilho" (cobertura do solo) e 
"verde" (crescimento da vegetacao), por moio da transf ormacao de 
"tasseled cap" (Kauth and Thomas, 1976); e a extracao de atributos. espa 
ciais por meio de filtros locais, que procuram levar em conta a informa 
cao da textura da imagem para a classificacao (Dutra e Mascarenhas, 
1980). 


Alen> disso o analista dispora, a exemplo do que acontece 
no LACIE (Spiers and Patterson, 1978), de informacao adicionala respej[ 
to dos pontos a serem rotulados, com a finaiidade de ressaltar carac 
rTsticas do desenvolvimento da cultura. Esta informacao consiste: (1) 

na aplicacao da transformacao de "tasseled cap" aos pontos de interes^ 
se e determinacao das coordenadas dc "brilho" e "verde"; (2)nadetermi^ 
nacao do espalhamento espectral desses pontos para cada passagem e da 
trajetoria espectral de cada ponto (valores de "brilho" e "verde" emca 
da passagem). 

2.2 - CLASSIFICACAO 

A operacao de classificacao separa cada ponto do segmento 
em duas classes: cultura "x" e "outros", de acordo com a regra de maxi 
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ma verosimilhanca. Adicionalmente, supoe-se que as classes cultura "x" 
e "outros" sao divididas em subclasses, espectrais, considerando-se que 
a distribuicao de cada uma delas seja gaussiana. 0 numero de subclasses 
e os pontos que compoem cada subclasse sao determinados atraves deumal 
goritimo de aglomeracao ("clust°ring") do tipo "vizinho-ma is - proximo" , 
que utiliza os pontos do tipo 1 como sementes. A parti r dos pontos de 
cada agregamento, sao calculados os parametros das distributes (media 
e matriz de covariancial de cada subclasse. As subclasses sao entao ro 
tuladas (cultura x e "outros") pelo criterio de minima distancia eucli_ 
diana. 0 rotulo da subclasse e o rotulo do ponto de tipo 1 mais proximo. 

Por fim, todos os pontos do segmento sao classificados se 
gundo a regra de maxima verossimilhanca . 

2.3 - ESTIMACAO DA PROPORCAO DE CULTURA NO SEGMENTO 

0 resultado da classificacao e usado para a estimacao de 
proporcao de cultura no segmento. Nesta fase, os "pontos-tcste" rotula 
dos pelo analista na fase de rotulacao-treinamento sao utilizados para 
corrigir eventuais tendencias Coias") na classificacao. 

De maneira simpl if icada , o procedimento de correcao de 
tendencia obtem o estimador p para a proporcao de cultura no segmento a 
partir da equacao: 

P = Pi: • * + Pio(1 ~ 0* (1) 

onde p n e a razao entre o numero de pontos do tipo 2, denominadcs clas 
se 1 (cultura "x") pelo analista e classificados como classe 1 pela ma 
quina, e o numero de pontos do tipo 2, classificados pelo computador co 
mo classe 1; e p 10 e a razao entre o numero de pontos do tipo 2, denonn 
nados classe 1 pelo analista e classificados como classe 0 ("outros") pe 
la maquina, e o numero de pontos do tipo 2, classificados como classe 0 pelo 
computador; \ e a razao entre o numero de pontos classificados como clas^ 
sificados como classe 1 pela maquina e o numero de pontos do segmento. 
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2.4 - AVALIACAO - APROVACAO 


0 resultado da estimegao e submetido ao anal ista para apro 
vag : io; a seu criterio, um ou mesmo todos os passos do procedimento pode 
rao ser rjfeitos. Para a avaliag^o sao tambem utilizados mapas de das 
sificagao. Uma vez obtida a proporgao de cultura para cacia segmento, os 
dados sao agregados para fornecer a area p 1s »ntada em uma determinada re 
giao. 

3. CONCLUSAO 


0 procedimento apresentado deve ser encarado como uma pH 
meira aproximagao de uma metodologia satisfatoria para estimativa de 
areas agrTcolas em imagens de satelite; com base na experiencia adquiH 
da, futuros desenvol vimentos procurarao melhorar o desempeniio do proce 
dirnento. Apontam-se aqui algumas areas critTcas, que certamente merece 
rao estudos posteriores: 

a) 0 aigoritmo de aglomeragao para obtengao as subclasses ("viz^ 
nho-mais-proximo") e muito simples para ser otimo. 0 desenvolH 
mento de algoritmos mais sofisticados, que utilizem informagao 
espacial juntamente :om a espectral , vievera melhorar o desempe 
nho da fase de classificagao. 

b) A rotulagao de pontos isolados pelo anal ista pode ser uma fonte 
de erro, devido a diversos fatores como a nistura de c. asses den 
tro de um ponto. Uma alternativa :eria fazer com que o anal ista 
rotulasse "regioes homogeneas" obtidas previamente atraves de 
aglomeragao. 

c) A hipotese gaussiana para a distribuigao das subclasses 20 algo 
ritmo de classificagao por maxima verossimilhanga podem nao ser 
os mais adequados para 0 probleina, dada a heterogeneidade das ca 
racterTsticas disponiveis. Algoritmos ma’S apropriados a nature 
za dos dados serao, certamente, objeto de futuras pesquisas. 
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